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Abstract 

It is well known that if a crystal structure has no inversion symmetry, it may allow for 
Dzyaloshinskii-Moriya magnetic interactions, operating between different crystallographic unit 
cells, which in turn should lead to the formation of long-periodic spin-spiral structures. Such 
a behavior is anticipated for two simple perovskites PbVOs and BiCoOs, crystallizing in the non- 
centrosymmetric tetragonal PAmm structure. Nevertheless, we argue that in reality PbVOa and 
BiCoOs should behave very differently. Due to the fundamental Kramers degeneracy for the odd- 
electron systems, PbVOa has no single-ion anisotropy. Therefore, the ground state of PbVOa 
will be indeed the spin spiral with the period of about one hundred unit cells. However, the even- 
electron BiCoOa has a large single-ion anisotropy, which locks this system in the collinear easy-axis 
C-type antiferromagnetic ground state. Our theoretical analysis is based on the low-energy model, 
derived from the first-principles electronic structure calculations. 
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I. INTRODUCTION 



Magnetic materials, crystallizing in the noncentrosymmetric structure, have attracted a 
great deal of attention. The lack of the inversion symmetry gives rise to the ferroelectric 
activity. If the latter property is combined with the magnetism, the system becomes mul- 
tiferroic, which has many merits for the next generation of electronic devices: for example, 
one can control the magnetization by applying the electric field and vice versa. The canon- 
ical example of such materials is BiFeOs, which possesses simultaneously high magnetic 
transition temperature (about 640 K) and high ferroelectric Curie temperature (about 1090 
K).l 

Recently fabricated PbVOs and BiCoOs belong to the same category. They crystallize 
in the noncentrosymmetric tetragonal PAmm structure (Fig. [T]).-i^ BiCoOs is an antiferro- 




FIG. 1. (Color online) Fragment of the crystal structure of BiCoOs. The Bi atoms are indicated 
by the big blue (dark) spheres, the Co atoms are indicated by the medium red (dark grey) spheres, 
and the oxygen atoms are indicated by the small green (light grey) spheres. 

magnetic (AFM) insulator of the C-type with the Neel temperature of about 470 K.^ The 
experimental information about PbVOa is rather controversial.-'- The magnetic susceptibil- 
ity has a broad maximum around 200, which might be the sign of the antiferromagnetism. 
On the other hand, no long-range magnetic order was found down to 1.8 K in the neutron 
diffraction experiments. However, the analysis of the neutron data depends on the model of 
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the magnetic structure, which is typically assumed in the process of interpretation. Finally, 
the experimental studies of PbVOa were hampered by possible defects in the sample.- 

According to first-principles electronic structure calculations, both PbVOs and BiCoOs 
are expected to have the C-type AFM ground state, although in PbVOa it is nearly degen- 
erate with the G-type AFM state.-^'^ Giant electric polarization (more than 150 /iC/cm^) 
was predicted theoretically both for PbVOs and BiCoOs,^ which spurred additional interest 
to these systems. 

Nevertheless, the violation of the inversion symmetry gives rise to a number of interesting 
effects, which are not currently accessible by the first-principles electronic structure calcula- 
tions, simply due to their complexity. One of them is a complex magnetic ordering, caused 
by antisymmetric Dzyaloshinskii-Moriya interactions: in the noncentrosymmetric systems, 
these interactions, of the relativistic origin, can operate between different crystallographic 
unit cells, thus driving the formation of long-periodic spin-spiral superstructures.- Particu- 
larly, the idea of the spin-spiral order in various oxide materials has attracted much attention 
recently in the context of their multiferroic behavior and was proposed as one of the possible 
origins of such behavior.- 

In this paper we will address some basic issues of the formation of the spin-spiral states 
in PbVOs and BiCoOs. We will argue that, despite similarities in the lattice distortion 
and population of the crystal-field levels, these two compounds will behave very differently. 
Particularly, we will show that the long-periodic spin spiral is a probable candidate for 
the magnetic ground state of PbVOs, where due to the fundamental Kramers degeneracy, 
the single-ion anisotropy does not exists. On the contrary, the spin-spiral state in BiCoOs 
(for which the Kramers theorem is no longer applicable) is suppressed by the single-ion 
anisotropy, which reinforces the formation of the easy-axis coUinear C-type AFM ground 
state. Our analysis is based on the low-energy model, derived from the first-principles 
electronic structure calculations. In this sense, this is the continuation of our previous works, 
devoted to 'realistic modeling' of complex oxide materials and other strongly correlated 
systems.—^— 

The rest of the paper is organized as follows. In Sec. Hi] we briefly discuss the construction 
of the low-energy model (in our case - the multiorbital Hubbard model) on the basis of 
first-principles electronic structure calculations. All model parameters can be found in the 
supplemental materials.— Sec. IIII Al is devoted to semi-quantitative analysis of the spin 
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model, which can be derived from the multiorbital Hubbard modeL Particularly, we consider 
the formation of incommensurate spin-spiral states, resulting from the competition of the 
isotropic exchange and Dzyaloshinskii-Moriya interactions, and explain the main difference 
in the behavior of the single- ion anisotropy in PbVOs and BiCoOa. In Sec. IIII B[ we will 
present results of extensive Hartree-Fock calculations for the long-periodic spin-spiral states 
in the electronic Hubbard model. Finally, in Sec. |IV] we will briefly summarize the main 
results. 



II. CONSTRUCTION OF THE LOW-ENERGY MODEL 



The magnetic properties of PbVOs and BiCoOs are mainly determined by the behavior 
of 3(i-bands located near the Fermi level. Therefore, our basic idea of our approach is to 
construct an effective low-energy model, formulated in the Wannier-basis for the Sd-bands, 
and to solve it by using model techniques. More specifically, we adopt the form of the 
multiorbital Hubbard model on the lattice of transition-metal sites: 

^ = XI + IY.Y1 Uaa'a"a"'cla4a"CiaAa"' ■ (l) 

ij aa' i {a} 

where we use the shorthand notations, according to which each Greek symbol stand for the 
combination of spin (s= t or i) a^nd orbital (m= xy, yz, 3z^— r^, zx, or x'^—y'^) indices. 
All parameters of the model Hamiltonian can be derived in an ah initio fashion, on the 
basis of first-principles electronic structure calculations. For instance, the one-electron part 
tf""' was obtained by using the downfolding procedure, and the Coulomb (and exchange) 
interactions Uaa'a"a'" ^ by combining the constrained density-functional theory (DFT) with 
the random-phase approximation (RPA). The method was discussed in the literature, and 



for details the reader is referred to Ref. 
can be found in Refs. 
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Resent applications to multiferroic compounds 
In all calculations we use experimental parameters of the 
crystal structure, reported in Refs. yandS 
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and 
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Without spin-orbit interaction, t^-^ is diagonal with respect to the spin indices t^" = 
t^"^'6ss'- The site-diagonal part of = \\t^^'\\ describes the crystal-field effects, while the 
off-diagonal part stands for transfer integrals. 

The crystal field stabilizes the xy orbitals (Fig. [2]). The splitting between xy- and the 
following after them yz- and zx-orbitals is about 1 eV, both for PbVOa and BiCoOs. The 
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FIG. 2. Scheme of the crystal-field splitting in PbVOs (left) and BiCoOs (right). 



x^—y^ and 82;^— orbitals lie in the higher-energy region (and are substantially higher for 
PbVOs, in comparison with BiCoOa). Thus, from the viewpoint of the crystal-field splitting, 
PbVOs and BiCoOs are expected to be very similar: PbVOs has only one ci-electron, which 
will occupy the xy orbital. Amount six d-electrons of BiCoOa, the one for the minority-spin 
states will also occupy the xy orbital and interact with the spherical ci-electron density of 
occupied majority-spin shell. 



The details of transfer integrals can be found in Ref. [iSl. One of the most interesting 
features in ti^ is the appearance of the so-called 'forbidden hoppings', for example between 
2^2 g^j^j orbitals in the bond 1-2, which would not exist in the centrosymmetric 
structure (see Fig. [1] for the notations of atomic sites). These transfer integrals have the 
following form (in meV): 
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for PbVOs and BiCoOs, respectively, in the basis of xy, yz, Sz'^—r'^, zx, and x'^—y'^ orbitals. 
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The microscopic origin of such forbidden hoppings was considered in Ref. 
parity violation, the Wannier orbital, which is formally labeled as "3z^— r^" has some weight 
of the Pz orbitals and, therefore, can interact with the zx orbitals of the neighboring sites. 
Such hoppings give rise to the antisymmetric part of tij, which is responsible for the ap- 
pearance of Dzyaloshinskii-Moriya interactions. Similar situation takes place in the bonds 
1-4 and 1-5. On the contrary, due to the rotational symmetry, the transfer integrals in the 
bond 1-3 are diagonal with respect to the orbital indices (the actual values are t^™= —45, 4, 
-237, 4, and -21 meV for PbVOg, and ^^^3™= -19, 24, -42, 24, and -34 meV for BiCoOg 
- other details can be found in Ref. Il3[ ). Therefore, the 'forbidden hoppings' do not take 
place and Dzyaloshinskii-Moriya interactions will vanish. 

For the relativistic spin-orbit interaction (SOI), we adopt two schemes. In the first one, 
we evaluate SOI only at the transition-metal sites and add it to the site-diagonal part of t^"' 
in the form ^LS, where ^=35 and 81 meV for PbVOa and BiCoOs, respectively. Thus, the 
effects of the SOI are expected to be larger in BiCoOs: due to large ^ and smaller crystal- 
field splitting (Fig. [2]), which competes with SOI. If it is not specified otherwise, we will refer 
to this scheme, for which most of the calculations have been performed. Nevertheless, as a 
test, we use also the second scheme, where SOI was included at all atoms on the level of the 
band-structure calculations and then corresponding parameters t^"' were derived through 
the downfolding procedure. For example, this schemes takes into account the effect of large 
SOI at heavy atoms Pb and Bi. 

The spin-dependence of Coulomb matrix elements has the standard form: Uaa'a"a"' = 



Umm'm"m"'^ss'^s"s"' ■ The details of f/, 
the strength of the matrix elements U, 



mm' m"m' 
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Rough idea about 



mm' m"m' 



can be found in Ref. 
can be obtained by interpolating them in 
terms of three characteristic averaged parameters f/, J and S, which would take place in 
the centrosymmetrical environment of isolated atoms. In these notations, U = F° is the 
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on-site Coulomb interaction, J = (F^ + F^)/14 is the intraatomic exchange interaction, and 
B = (9F^ — 5F^)/441 is the 'nonsphericity', in terms of radial Slater's integrals 
and F^. In the other words, U enforce the charge stability of certain atomic configurations, 
while J and B are responsible for the Hund rules. The results of such interpolation are 
shown in Table [H One can clearly see that the on-site Coulomb repulsion U is strongly 



TABLE I. Averaged values of the Coulomb interaction U, exchange interaction J, and the non- 
sphericity B, obtained from the fitting of the matrix elements Umm'm"m"' ■ All parameters are 
measured in electron volt. 



compound 


U 


J 


B 


PbVOa 


1.57 


0.84 


0.08 


BiCoOg 


2.38 


0.90 


0.09 



screened, especially in PbVOs, while other parameters are close to atomic values. We use this 
interpolation only for explanatory purposes, while all practical calculations were performed 
with actual parameters Umm'm"m"' reported in Ref. Il3|. The deviation of Umm'm"m"' from 
the centrosymmetric form is quite strong. For example, in the case of BiCoOa, the diagonal 
matrix elements vary as Ummmm= 3.84, 3.39, 2.94, 3.39, and 3.48 eV for m= xy, yz, Sz'^—r'^, 
zx, and x'^—y'^, respectively. 



After the construction, the model is solved in the Hartree-Fock (HF) approximation.- 
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III. RESULTS AND DISCUSSIONS 



A. Qualitative analysis based on the spin Hamiltonian 

The existence of the spin-spiral states in noncentrosymmetric perovskites can be under- 
stood in the framework of the spin model:— 

(ij) (ij) * 

(where Jij is the isotropic exchange interaction, dij is the vector of antisymmetric Dzyaloshinskii- 
Moriya interactions, and fa is the single- ion anisotropy tensor), which can be obtained by 
mapping the electronic model ([1]) onto the spin one and integrating out all degrees of free- 
doms but spins. There are several ways how to do it. One possibility is to consider the 
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perturbation-theory expansion with respect to the infinitesimal spin rotations and SOI near 
the nonrelativistic ground state in the Hartree-Fock approximation.— In the following, the 
results of such model will be denoted by the symbols 'm/'. Moreover, for the d} config- 
uration of PbVOs one can easily consider the theory of superexchange interactions in the 
second order with respect to the transfer integrals (in the following denoted by the symbol 
'set').- 

Then, neglecting for a while the single-ion anisotropy term, the energy of (classical) spin 
spiral in the zx-plane, 

{Si) = 5 (sinqRj, 0, cosqRj) 
(Rj being the radius- vector of the site i), is given by 

E{q) = - ^ ( Joi cos qRi - dl^ sin qRj) , 

i 

and the spin-spiral vector q = (gx, vr, 0) in the ground state should correspond to the min- 
imum of E{q). Obviously, the isotropic exchange interactions (Joi) will tend to establish 
a collinear spin structure with qRj= or tt, while Dzyaloshinskii-Moriya interactions {d^j) 
will deform this structure and make it incommensurate.— The Dzyaloshinskii-Moriya interac- 
tions are also responsible for the asymmetry between right-handed {q^ > 0) and left-handed 
{qx < 0) spin-spiral states, which is manifested in the inequality -E'(q) 7^ E{—q). 

Parameters of isotropic exchange interactions {Jij) are listed in Table [Tll We note that 

TABLE II. Isotropic Heisenberg interactions (measured in meV) for PbVOs and BiCoOs. Notations 
of the atomic sites are explained in Fig. [TJ Results of the superexchange theory are denoted by 
the symbols 'set'. Results for infinitesimal spin rotations near the nonrelativistic ground state are 
denoted by the symbols 'in/'. 



bond PbVOs (set) PbVOs (m/) BiCoOg (inf) 

1-2 -49.86 -44.71 -9.65 

1-3 -3.63 -0.64 -0.15 

1-4 4.76 3.20 -0.91 

1-5 3.94 1.25 -1.19 

1-6 3.61 1.25 -0.07 



the schemes 'set' and 'm/' in the case of PbVOa provide very similar results. This seems to 
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be reasonable, because if the orbital configuration is quenched by the crystal-field splitting, 
the spin model ([2]) is well defined and the parameters are not sensitive to the way how they 
are defined (of course, provided that \tij/U\ « 1 and the schemes 'set' makes a sense). The 
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magnetic transition temperature, evaluated in the random-phase approximation (see Ref. 
for details) for the G- and C-type AFM states, is of the order of 200 and 600 K for PbVOs 
and BiCoOs, respectively. The experimental Neel temperature for BiCoOs is 470 K.^ The 
situation in PbVOs is rather controversial. On the one hand, the results of the neutron pow- 
der diffraction experiment are not conclusive, because their interpretation strongly depends 
on the magnetic structure, which was assumed for the analysis of experimental data.^ On 
the other hand, the magnetic susceptibility of PbVOa does display a broad maximum at 
around 200 K, which could be regarded as the sign of an antiferromagnetism. Moreover, the 
G-type antiferromagnetic order was proposed for the thin films of PbVOs below 130 K.— 
Parameters of Dzyaloshinskii-Moriya interactions are shown in Table lllli They are at 

TABLE III. Nonvanishing parameters of Dzyaloshinskii-Moriya interactions djj = (df,-, dfj, df^) 
(measured in meV) for PbVOs and BiCoOs. Notations of atomic sites are explained in Fig. [H 
Other parameters are equal to zero. Results of the superexchange model are denoted by the 
symbols 'set'. Results for infinitesimal spin rotations near the nonrelativistic ground state are 
denoted by the symbols 'm/'. Note that Dzyaloshinskii-Moriya interactions vanish in the bond 
1-3 due to symmetry constraints. 

parameters PbVOs (set) PbVOs {inf) BiCoOs (inf) 

42 -0.98 -0.77 -0.13 

dl^ = -dli 0.39 0.17 

^5 -0.11 -0.04 

4g = -dfg -0.01 -0.03 



least one order of magnitude smaller than for the same bonds. 

Using these parameters, the spin-spiral vector q^, can be estimated as g^^a = tt — A<^ 
(a being the lattice parameter in the xy-plane), where A0= 6 x lO^vr and 4 x IO'^tt for 
PbVOs and BiCoOs, respectively. Thus, by considering only Jij and djj, both materials are 
expected to form spin-spiral structures, involving more than one hundred unit cells. As we 
will see below, this scenario indeed holds for PbVOs, but not for BiCoOs. 
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The main difference between PbVOs and BiCoOs is in the behavior of the single-ion 
anisotropy fa. For the S = 1/2 compound PbVOs, hi is expected to be zero as the con- 
sequence of fundamental Kramers degeneracy for the odd-electron systems. Particularly, 
the ground state of the self- interaction free ion V^^ is the Kramers doublet. Therefore, 
the rotation of spin corresponds to the unitary transformation of the wave function within 
this doublet without any energy cost. The situation is completely different for the S = 2 
(or even-electron) compound BiCoOs: the Kramers theorem is no longer valid, which for- 
mally allows for the finite fa. This statement can be verified by direct calculations of the 
anisotropy energies AE = — Ej_ (where the symbols "||" and "_L" correspond to the spin 
configurations, where (Si) is parallel and perpendicular to the tetragonal z-axis). In the 
C-type AFM state, it yields AE= 0.02 and —5.63 meV per formula unit for PbVOs and 
BiCoOs, respectively. Moreover, the main contribution to AE indeed originates from the 
single-ion anisotropy. This can be seen by repeating the same calculations in the atomic 
limit (and enforcing tij=0 for all i^j), which yields AE= and —5.86 meV per formula 
unit for PbVOs and BiCoOs, respectively. Small deviations from the atomic limit are due to 
intersite (i^j) anisotropic interactions fij, which can be evaluated in the 'set'-model and are 
at least one order of magnitude smaller than dij}^ Using the obtained values of AE and the 
symmetry considerations, nonvanishing parameters of the single-ion anisotropy for BiCoOs 
can be estimated as t^^ = t^^ = —\t^^ = 0.49 meV. Thus, we are dealing the following 
hierarchy of magnetic interactions | Jij|>>|fjj|>>|djj|. It means that the formation of the 
spin-spiral state in BiCoOs is strongly affected by the single-ion anisotropics, which will 
tend to restore the collinear spin structure by aligning the magnetic moments either parallel 
or antiparallel to the 2;-axis. Of course, the final answer about the from of the magnetic 
ground state of BiCoOs can be obtained only on the basis of detailed calculations, which we 
will discuss in the next section. 

B. Solution of electronic model 

In this section we present results of extensive Hartree-Fock calculations for large super- 
cells, which allow for the spin-spiral solutions with q^a = 7r(|L| — 1)/L, where \L\ is the 
number of cells along the x-axis: L > and < corresponds to the right- and left-handed 
alignment, respectively, and the limit |L| — )■ oo corresponds to the collinear C-type AFM 
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state. The main results are summarized in Figs. [3] and |H One can clearly see that there 



^ ^ ^ (a) 



\ \ I i / ^ ^ \ \ I ^(b) 



FIG. 3. Spin magnetic structure along the x-directions in the case of PbVOa (a) and BiCoOs (b), 
as obtained in the Hartree-Fock calculations for L= 21. Here, x is the horizontal axis and z is the 
vertical one. 





FIG. 4. Distribution of the spin magnetic moments (as if they were brought to the same origin) in 
the xz-plane of PbVOs (top) and BiCoOs (bottom): results of Hartree-Fock calculations for L= 
7 ('a' and 'e'), 21 ('b' and 'f'), 35 ('c' and 'g'), and 49 ('d' and 'h'). Here, x is the horizontal axis 
and z is the vertical one. 



is a big difference between PbVOs and BiCoOa. PbVOs tends to form a homogeneous 
spin-spiral state, where the angle between neighboring magnetic moments along the x-axis 
remains constant (small deviations are caused by weak inter-site anisotropy effects). On 
the contrary, due to the large single-ion anisotropy, the spin-spiral configurations in BiCoOs 
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are strongly distorted, and the moments are bunched around the z-axis (Figs. Hf-h). The 
so-called 'bunching effect' is well know for magnetic rare-earth metals and was intensively 
discussed already more than forty earth ago.— Thus, BiCoOa tends to form an inhomoge- 
neous magnetic state, which corresponds to the (nearly) collinear AFM alignment in the 
wide part of the supercell, except small 'domain wall', where the spins undergo the reorien- 
tation within the area of about ten unit cells. The latter solutions were obtained for odd 
numbers of cells L, which in the AFM lattice results in the formation of the domain wall 
defect. For even L, the Hartree-Fock equations converge to the C-type AFM state. The spin 
pattern in the domain wall is well reproduced already for L = 21 (Fig. Hj). For larger cells, 
the additional spins participate in the formation of the AFM regions, and are either parallel 
or antiparallel to the 2;-axis, leading to tiny changes in Figs. Hf-h, which are practically not 
distinguishable to the eye. 

Results of total energy calculations (Fig. [5]) are well consistent with the above finding. 
As expected for the spin-spiral states, the dependence of the total energy on 1/L in the case 
of PbVOs is well described by the parabola. Due to the Dzyaloshinskii-Moriya interactions, 
there is a small asymmetry of the total energy with respect to the inversion L — )■ —L of 
chirality of the spin spiral. Thus, the total energy minimum, obtained from the extrapola- 
tion, corresponds to the spin-spiral ground state with L ^ 94. On the contrary, the total 
energy of BiCoOs is a linear function of This is because of the localized character of the 
domain wall, for which the total energy (divided by the total number of cells) is expected 
to scale as 1/L. Thus, the minimum corresponds to the collinear C-type AFM ground state 
(|L| — i- oo), in which the total energy exhibits the derivative discontinuity. Nevertheless, 
even in this case, the total energy has different slops in the regions L > and L < 0, again, 
due to the difference between the right- and left-handed spin-spiral alignment in the domain 
wall. 

The crucial role of the single-ion anisotropy in the formation of the easy-axis C-type AFM 
ground state can be illustrated by repeating supercell calculations for BiCoOa with the same 
parameters of the model Hamiltonian, but with different number of valence electrons: one 
instead of six. Thus, according to the Kramers theorem, the single-ion anisotropy should 
vanish, similar to PbVOs. The results for L = 21 are shown in Fig. [6], in comparison 
with regular BiCoOs, including all six valence electrons. There are two effects. First, 
as was already discussed in Sec. [Tll the effects of SOI are generally larger in BiCoOs, in 
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FIG. 5. (Color online) Total energies versus 1/L as obtained in the Hartree-Fock approximation 
for PbVOs (left) and BiCoOs (right). Calculated points are denoted by symbols. Solid line is the 
results of interpolation E = a/ L + h/ L? in the case of PbVOs (where a = —0.99 meV and h = 46.36 
meV) and E = a/L in the case of BiCoOs (where a = —42.22 meV for L < and a = 36.94 meV 
for L > 0). The location of the total energy minimum is shown by the dot-dashed line. 

comparison with PbVOa. Therefore, inter-site anisotropic interactions become stronger, 
which is reflected in some bunching of the spin magnetic moments around the horizontal 
X-axis in the hypothetical 'single-electron BiCoOa' (similar bunching exists in PbVOs - 
Fig. m but the effect is considerably weaker). Second, the easy-axis alignment in BiCoOs is 
solely related to the single-ion anisotropy term: as long as it is absent in the hypothetical 
'single-electron BiCoOs', the spin magnetic moments start to regroup around the x-axis. 

Finally, we comment on the dependence of our results on different levels of treatment of 
the relativistic SOI. We consider two such schemes: (i) the SOI was included to the model 
Hamiltonian as a pseudo-perturbation only at the transition-metal sites," and (ii) the SOI 
was included at all sites of the system (including heavy Pb and Bi elements) in the process of 
downfolding procedure. However, the distribution of the spin magnetic moments, obtained 
in these two schemes, is practically indistinguishable (Fig. [7]). Thus, the SOI at the heavy 
Pb- and Bi-elements does not seem to play an important role in the magnetic properties of 
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FIG. 6. Distribution of the spin magnetic moments (as if they were brought to the same origin), 
obtained in the Hartree-Fock calculations for L= 21. Left panel (a) shows results for the regular 
BiCoOs, involving six valence electrons. Right panel (b) shows the same results for the hypothetical 
system, which has the same parameters of electronic Hamiltonian as for BiCoOs and only one 
valence electron. Thus, according to the Kramers theorem, the single-ion anisotropy terms should 
not operate in the case (b). Here, x is the horizontal axis and z is the vertical one. 




FIG. 7. Influence of the spin-orbit interaction on the magnetic structure of BiCoOa (results for 
L= 21). Left panel (a) shows the results, where the spin-orbit interaction was included only on 
the Co-sites. Right panel (b) shows the results, where the spin-orbit interaction was included on 
all sites of the system, through the downfolding procedure. Here, x is the horizontal axis and z is 
the vertical one. 



PbV03 and BiCoOs. 
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IV. CONCLUSIONS 



Being based on results of the low-energy electronic model, derived from the first-principles 
electronic structure calculations, we analyzed possible magnetic structures of two noncen- 
trosymmetric perovskites PbVOa and BiCoOa. We have argued that, despite structural 
similarities, the magnetic behavior of these two materials is expected to be very different. 
PbVOs, with the spin S = 1/2, should form a long-periodic spin-spiral state, which results 
solely from the competition between isotropic exchange and Dzyaloshinskii-Moriya inter- 
actions in the noncentrosymmetric crystal structure. Due to the Kramers degeneracy, the 
single- ion anisotropy does not operated in PbVOa. However, the latter is expected to play 
a major role in BiCoOs, which has the spin S = 2. Particularly, the single-ion anisotropy 
suppress the noncoUinear spin-spiral alignment in BiCoOs and enforces the formation of the 
C-type antiferromagnetic ground state, in agreement with the experiment.™ 

We believe that this funding has a direct implication to the properties multiferroic man- 
ganites, which also have spin S = 2 and the large single- ion anisotropy. Therefore, the 
numerous claims about the spin-spiral ground state of these compounds, and related to it 
improper ferroelectric activity, should be taken cautiously. Again, due to the large single-ion 
anisotropy, the ground state of manganites is not necessary the spin spiral, which prompts 
a search for alternative mechanism of multiferroicity in these compounds.— 
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